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INTRODUCTION 


The atomic iodine laser based on the photodissociation of alkyl iodide gases was 
among the first proposed for direct solar- to-laser energy conversion in space 
(ref. 1). The first solar-simulator demonstration being achieved (ref. 2), more 
attention is now given to system requirements of chemical reversibility and lasant 
reprocessing in scaling to large space systems. Such a scaling requires an adequate 
knowledge of the chemical processes present in the gas during laser operation 
(ref. 3). The development of a preliminary chemical kinetic model is presented 
elsewhere with limited comparison to solar simulator experiments (ref. 4) . 

It was shown in a previous report (ref. 4) that the laser output consists of a 
transient pulse region followed by a continuous wave (cw) region in which the laser 
power output is proportional to the power absorbed. At late times, the cw region is 
terminated as a result of quenching by the iodine molecule which formed through 
three-body collision of two ground-state iodine atoms with an alkyl iodide molecule. 
The I 2 formation critically depends on the reaction rates of the alkyl iodide 
recombination 


R + I ■> RI 


and three-body formation of I 2 as 


I + I + RI -*■ I 2 + RI 


where R denotes the alkyl radical. It was determined (ref. 4) that some of the rate 
constants for these reactions in the literature are not compatible with Langley 
experiments. 

In the present report, the kinetic factors involved in reaching laser threshold 
are examined. In particular, the effects of gas diffusion, hyperfine splitting, and 
gain length scaling are treated in detail. 


SOLAR SIMULATOR 

The solar-simulator light source is an arc discharge across an 8-mm gap stabi- 
lized by a 1 030-kPa (10.2-atm) Xe flow. The light output is reflected by a high- 
quality elliptical reflector and focused in the chopper plane as shown in figure 1. 
The reflector has a vapor-deposited aluminum front surface protected by a MgF 2 coat- 
ing to prevent surface oxidation and abrasion. The spectral content of the arc 
plasma output corresponds to blackbody emission at 6000 K at the focus of the ellip- 
tical reflector. The divergent light beam leaving the focus is collected by a 
7075 aluminum alloy 25° cone with a polished interior surface and MgF 2 coating. The 
light beam is again focused on the cone axis, which is aligned with the optical axis 
of the system. Three different cone sizes were used to gain scaling information for 
model testing. 



50 kW 

Solar Simulator 



Figure 1.- Setup for solar-simulator-pumped laser experiment. 


Hie photon flux on the wavelength interval between A and A + dX emitted by a 
blackbody radiator at a temperature T is 


<t>( X) d A. = 


2tec dX 

4 exp(hc/XkT) 
X 


1 


(1 ) 


where c is the velocity of light, h is Planck's constant, and k is Boltzmann's 
constant. The exponential argument is given as 


he 1 .44 
kAT “ XT 


( 2 ) 


for X in centimeters and T in kelvin. Hie light is reflected from the conic 
collector with a reflectance R^ at each wavelength before reaching the optical 
axis. Hie light intensity at the optical axis is related to the geometry of the arc 
plasma and the reflector and collector surfaces, denoted by the function M(r, Jt), as 


I (r,Jl) dX = M(r,i) R e <j> ( X) dX 
X XX 


(3) 


where e is the photon energy. 

X 

Hie light intensity is measured by a calorimeter and represented by 


C (r , JO = exp 


U 2 r 2 ' 

■ 2 - 77 N 

\L Rl 


(4) 


2 



where CM is the peak intensity of the arc image, L is the full length at half 
maximum light intensity along the tube axis, R is the full width at half maximum 
light intensity transverse to the tube axis, A is the distance along the tube axis, 
and r is the distance from the tube axis. The size of the arc image depends on the 
size of the cone collector and its distance from the focus. It is clear that 


C(r , A) 



M(r ,X) R £^ 


<t>(\) dX 


(5) 


which allows M(r,A) to be determined from the calorimeter data. Since most of the 
energy is in the visible region where R^ is nearly constant, equation (5) may be 
written with the help of equation (4) as 


M (r , A) 


CM exp [-2.7 7 (A 2 /L 2 + r 2 /R 2 )] 


R ( £ <t>(\) d\ 

vis Mq a. 


(6) 


where is replaced by the constant R v i s , since mainly visible radiation contrib- 

utes to the calorimeter measurements. Expressing the numerator and denominator in 
units of the solar constant yields 


M(r , A) 


C q exp [-2.7 7 (a 2 /L 2 + r 2 /R 2 )] 


R 

vis 


(7) 


where it is found from the calorimeter data that 


2.72 x 10 4 

2 

(1 + 0.375m) 


( 8 ) 


L 


( 4 + 1 . 5m) (m + 5) 
8 


(9) 


2.19 2 + 0.824m 
5 + m 


( 10 ) 


and m is a parameter associated with the scale of the three collector cones with 
different apertures. (See table I.) 
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TABLE I 


SCALE PARAMETERS FOR COLLECTOR CONES 


Aperture 

diameter, 

cm 

Distance from 
focus, cm 

m 

C o' 

solar constants 

L, 

cm 

R» 

mm 

9.6 

7.0 

0.5 

19 288 

3.27 

4.73 

16.6 

15.2 

2.0 

8 881 

6.13 

5.49 

23.6 

24.5 

2.5 

7 245 

7.27 

5.67 


A comparison between measurements of intensity along the optical axis and the 
model function of equation (4) in units of the solar constant is shown in figure 2. 
Measurements in the radial direction for the cone with an aperture diameter of 
16.6 cm at a location 3 cm along the optical axis from the peak intensity are shown 
in figure 3 in comparison with a theoretical curve based on equation (4). 



Figure 2.- Illumination on axis of laser tube as a function of distance 
behind cone face. Curve is based on equation (4) and table I. Dots 
are e>q>erimental measurements of April 13, 1981 (aperture with 9.6-cm 

diameter). Triangles are experimental measurements of July 30, 1981 
(aperture with 16.6-cm diameter). Squares are experimental measurements 
of March 25, 1981 (aperture with 23.6-cm diameter). 
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• Experiment 
— Theory 


5000 1 



r, cm 

Figure 3.- Illumination as a function of distance from cone axis 
(1 6.6-cm-diameter aperture). Curve is based on equation (4) 
and table I. Dots are experimental measurements of July 30, 1981. 


Although the lamp intensity is presented as a static quantity, there are, in 
fact, a number of poorly understood factors that affect the illumination of the cen- 
tral region of the collector cone. Because the lamp electrodes are mechanically held 
in place, slow erosion will widen the gap and cause a less localized arc. Ihe 
installation of new lamps will not restore the machine to its initial operation 
because of lamp differences. In addition, there are unknown factors causing day-to- 
day variations in the measured illumination. Aside from these long-term variations, 
there are small temporal fluctuations in the arc plasma which cause large spatial 
movements in the arc image because of the large magnification factors. Ihis produces 
some difficulty in interpretation of the laser experiments , since the instantaneous 
pump power is generally less than its maximum. Only by considering many experimental 
runs can the ideal experimental result (that is, result corresponding to maximum pump 
power) be identified. 

The light is reflected by the conic collector and concentrated on the cone axis 
with intensity given by equation (3). The laser tube is centered on the cone axis 
for which the peak light intensity crosses the tube wall at an angle of 46° before 
reaching the lasant gas. (See fig. 4.) Ihe reflection loss at the tube interfaces 
is R » 0.08, and the transmission through the 1-mm quartz tube is near unity. Ihe 
attenuation of the light in passing through the filling gas to the laser tube center- 
line is 


<6 (x) = exp[-n cr(\) x] R M(0,1) R <t>(\) d\ (11) 

A. t A. 
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r id^ 



x ~ r id /cos 46 ° 


T 


laser tube 

Figure 4.- Typical light ray reflected into laser tube. 


where x is the slant distance (see fig. 
surface along the light path (x » r id / cos 
specific absorbing gas molecule, and 


4) from the centerline to the inner tube 
46°), n is the number density of the 


ff(X) 


a 

o 


exp 


-(x - x y 

o 


45 


(1 2 ) 


is the specific photoabsorption cross section with peak cross section a , line cen- 
ter X Q , and line width 5 Q given in table II according to references 5 to 7. Of 
the photons absorbed in the spectral region of equation (12), the absorbing molecules 
are fragmented to yield atomic iodine. The two alkyl iodides and I 2 photodissociate 
into fractions 3> of metastable iodine with values of 100 percent and 51 percent, 
respectively. (See table II. ) 

TABLE II.- PHOTOABSORPTION PARAMETERS USED FOR THE EQUIVALENT POWER OF 
ONE SOLAR CONSTANT EXPOSURE, 1.4 kW/m 2 


Parame ter 

n-C 3 F 7 I 

i-C 3 F 7 I 

*2 


(ref. 5) 

(ref. 6) 

(ref. 7) 

a , cm 2 

VO 

X 

o 

1 

VO 

6.2 x 10 -19 

VO 

• 

— 1 

X 

o 

1 

vo 

o 




X , nm 

272 

275 

499 

o 




5 , nm 

12.7 

14.5 

23.0 

o 




i> 

I* 

1.0 

1.0 

0.51 

S, sec -1 

3.04 x 10-3 

3.37 x TO " 3 

3.38 x 10- 2 

R X / R vis 

0.8 

0.8 

1 .0 

o 





f 


0.652 


0.653 


0.673 




The rate at which 


I* is formed in the laser tube is 


f 

£(r,£) = I n cr(X) R M(r,JL) R. expt-n a(X) x] <|>(X) dX 

I X* t A. 

•'n 


= ® IJt nR 


t M(r, Ji.) I 

•'n 


R a(X) e>q)[-n a(X) x] <j>(X) dX 
A. 


(1 3) 


which is evaluated through approximation with Gauss-Hermi te quadrature (ref. 8) and 
yields 


S(r,A) 


26 $ no R R. M(r,Jl) 
o I* o t X 

o 


{ 1 . 1 8 <)> (X ) 
o 


exp(-nc x) 
o 


+ 2 .9 5 [ 4>( X + 2.456 ) + (|> ( X - 2.456 )] exp(-0.223ncr x)} (14) 

o o o o o 

Equation (14) may now be simplified as 


R R, 

t X 

^(r,A) = — Sn[f e3q3(-no x) + (1 

R . O 

VIS 


f) exp ( -0. 2 23na Q x) ] C(r,Z) 


(1 5) 


where S (maximum photodissociation rate) and f (fractional absorption near line 
center) are given in table II. Hie usefulness of the Hermite quadrature is demon- 
strated in figure 5 where equation (15) is compared with an accurate numerical evalu- 
ation of equation (13) and the usual exponential approximation. 


ENERGY LEVELS AND REACTION PRODUCTS 

There are several chemical /physical processes important to the formation and 
loss of excited iodine. Several requirements must be met in order that specific 
reactions may occur. The first such requirement is that the change in internal 
energy be less than or on the order of the average thermal energy. 
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Figure 5.- Relative attenuation of light as a function of penetration 

into lasant gas (n-C-jF^U* 


Molecular Fragmentation 

The approximate bond energies of the alkyl iodide constituents are given in 
table III. These energies are mean molecular values compiled by Huheey (ref. 9) . 


TABLE III.- SINGLE- AND DOUBLE-BOND ENERGIES AND DISTANCES 
OF ALKYL IODIDE CONSTITUENTS 

[From ref. 9] 


Bond 

Bond energy, 

D 0 ' eV 

C-C 

3.58 

C=C 

6.24 

F-F 

1.60 

I-I 

1 .54 

C-F 

5.03 

C-I 

2.21 

I-F 

2.88 


Bond distance, 
r, pm 

154 

1 34 

142 

267 

1 35 

21 4 

191 
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First consider the ways in which the basic CF 3 CF 2 CF 2 I molecule could fragment and the 
corresponding enthalpies as shown in table IV. The fragments are ordered according 


TABLE IV.- CF 3 CF 2 CF 2 I FRAGMENTATION PRODUCTS AND 
CORRESPONDING ENTHALPIES 


Reaction product 

AH, eV/mol 

cf 4 + cf 2 cfi 

0.92 

C 2 f 4 + CF 3 r 

.92 

CF 2 CFCF 3 + IF 

1 .70 

•CF CF CF + I 
2 2 3 

2. 21 

•CF + *CF CF I 
3 2 2 

3. 58 

•CF CF CF I + F 
2 2 2 

5.03 

cf 2 cfcf 2 i + f 2 

5.80 


to increasing enthalpies so that the most likely fragments appear at the top of the 
table. It is clear that none of these products should be present during ordinary 
room temperature chemical processes, although photofragmentation or other energetic 
processes could potentially form some of the listed fragments. There are no absorp- 
tion bands in the visible region (ref. 10) which could give rise to the products 
listed. Thus, they could be directly formed only at pump rates with heat release 
sufficient for pyrolytic decomposition. Many of these fragments have been observed 
in pyrolysis studies (ref. 11). It is anticipated that pyrolysis will not be 
important in the present study, since relatively small energy deposits are sufficient 
to excite the gas to threshold. Photofragmentation will, of course, be of major 
importance. 


Photodissociation Products 

Ihe lowest order photodissociation process in the perf luoroalkyl iodides corre- 
sponds to the electronic excitation of a nonbonding electron of the iodine to an 
antibonding (it*) orbital (ref. 10). A similar transition occurs in the fluorine but 
only for hard ultraviolet (UV) photons. Ihe potential energy curves along the C-I 
bond direction are shown in figure 6. The curves were drawn using the Morse function 
for the ground-state potential with the spring constant according to Badger's rela- 
tion (ref. 12). The repulsive antibonding potential is approximated by 

V (r) = e* + D + A exp[-a(r - r ) ] (16) 

no c 
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> 

o 

> 


Figure 6.- Potential along R-I bond direction. 

where e* = 0.94 eV is the l(^P-|/2) excitation energy; D Q is the C-I bond energy; 
r is the C-I separation; and the parameters A, a, and r c are determined from the 
photodissociation parameters X and 6 in table II. The photodissociation is 
predominantly into the iodine metastable state (I*) (ref. 13) 

RI + h v -> R + I* 
o 

where v Q is the central frequency, and the potential energy between the fragments 
immediately after photoabsorption is »1 .39 eV. The resultant kinetic energies of the 
R and I* fragments are 0.60 and 0.79 eV, respectively. It is not anticipated that 
the photodissociation of I 2 will ever play a major role in production of I or I*, 
since I 2 is a serious quenching molecule and quickly halts laser operation. Clearly, 
the major constituents in addition to the lasant gas to be found at any time in the 
laser tube will be the photolysis products and subsequent reaction products. The 
possible reactions of the photolysis products must now be examined. 

Primary Reactions 

The possible reactions of the lasant molecules and the photodissociation prod- 
ucts are given in table V in the order of increasing enthalpies. The last two reac- 
tions of table V are unimportant because of the large energy input required, and the 


I* + I* •> I 

2 

reaction is omitted because recombination is along a highly repulsive potential sur- 
face. The exciplex formation 


I* + RI ■> RI* 
2 

10 



r, nm 




TABLE V.- PRINCIPAL REACTIONS OF PHOTODISSOCIATION 
PRODUCTS AND CORRESPONDING ENTHALPIES 

[R = *C F ] 

3 7 


Reaction 

AH, eV/mol 

R + R + R 2 

-3.58 

1* + 1* > i 2 

-3.42 

1* + 1 -> i 2 

-2. 48 

R + I a- RI 

■ -2.21 

I + ! > ! 2 

-1.54 

R + RI a- R 2 + I 

-1.37 

R + I* -> RI 

-1 . 27 

I* * I 

-.94 

RI + RI •> R 2 + l 2 

-.70 

R + RI -> R 2 4* I* 

-.43 

I* + RI ■> RI* 
2 

-.07 

I 4- RI -> R 4* I 2 

.67 

I* + RI ■> R + I 2 

1 .61 


is contained in the usual quenching coefficient discussed subsequently. The remain- 
ing reactions are all exothermic, and their corresponding rates are mainly limited by 
steric factors and activation thresholds. For example, the reaction among lasant 
molecules 


RI + RI ■> R + I 
2 2 


is first highly limited by the complicated structure of the RI molecule but, in addi- 
tion, has a threshold approaching 4 eV. In fact, this reaction seems to play, at 
most, a minor role in pyrolytic decomposition (ref. 11). 

Hie recombination 


R + R ■* R 

2 
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has no threshold and is mainly delayed by steric factors, as is 


R + I > RI 


The excess kinetic energy in both reactions is absorbed by the internal energy of the 
alkyl radicals so that a third body is not required. 

The important reaction for lasers 


R + I* -> RI 


is inhibited by the repulsive potential (fig. 6) between them. This reaction is the 
principal loss of excited iodine and the main limit in reaching laser threshold. 

Loss of the electronic excitation energy through various collisional processes 
denoted by 


I* -> I 


is the ultimate limiting factor in the laser operation. Such losses occur through 
collisional energy transfer, exciplex formation, and collision with the wall. Radia- 
tive decay plays a minor role because of the long lifetime of this metastable state. 

The formation of molecular iodine is a three-body reaction; the third body is 
needed to stabilize the reaction by absorbing kinetic energy to extract the excess 
energy of translational motion. These reactions 


I + I ->■ I 

2 


and 


I* + I * I 

2 


are particularly important to laser operation because the high quenching rate of the 
iodine molecule is the main limit to laser pulse duration and processing of used 
lasant material. 
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Secondary Reactions 


The possible secondary reactions of the new species formed through the primary 
reactions are shown in table VI. Of these, only two reactions appear to be of any 
potential importance. 'Hie first, 


R + I 


2 


•* RI + I 


is suppressed by steric factors as well as activation energy on the order of 1 eV. 

If this reaction occurs, it is probably only among hot radicals produced in the pho- 
todissociation of the lasant molecules. The second reaction. 


I* + R -*■ RI + R 
2 


is inhibited by an approximate 2-eV threshold and poor electronic- to-vibration cou- 
pling. This reaction and subsequent reactions in table VI are unimportant. 


TABLE VI.- SECONDARY REACTIONS AND 
CORRESPONDING ENTHALPIES 


[R = *C 3 F 7 ] 



Reaction 

AH, eV/mol 


R + I 2 a- RI + I 

1 

o 

• 


I* + R 2 -> RI + R 

.43 


^2 ^2 2RI 

.70 


I + R 2 RI + R 

1.37 


CHEMICAL REACTION RATES 

The specific reaction rates of the reactions which appear possible under simple 
energy considerations must be known before a complete chemical kinetic model can be 
derived. Specific reaction rate data are collected and summarized in this section. 
The symbols for the reaction rate coefficients are defined in table VII. 
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TABLE VII.- MEAN REACTION RATE COEFFICIENTS AND ASSOCIATED UNCERTAINTY FACTORS 

BASED ON LITERATURE VALUES 

The factor in parentheses gives the uncertainty limits associated! 

with the coefficient I 


Reactants 

Products 

Reaction rate coefficient, (cm 3 ) n /sec 






Symbol 

R = n-C 3 F 7 

R = i-C 3 F ? 

R + R 

r 2 

k 3 

2.6 x 1 0 _1 2 ( 4 ) - 1 

9.0 x 1 0* 1 3 (3.8)* 1 

R + I 

RI 

K 2 

2.3 x 10" 1 1 ( 3. 5) - 1 

3.9 x 1 0* 1 1 (4.3) ±1 

I + I + RI 

I 2 + RI 

C 2 

8.5 x IQ" 32 (S.S)* 1 

8.3 x io* 32 (S.S)* 1 

I + I + I, 

1 2 + I 2 

C 4 

3.8 x TO -30 ( 1 . 3) - 1 

3.8 x io* 30 (1 .3)* 1 

R + RI 

r 2 + I 

K 4 

3 x io‘ 16 

3 x 1 0* 1 6 

R + I* 

RI 

K 1 

5.6 x io -13 (6 . 2) - 1 

1.7 x 1 o* 1 3 (1 7 ) “ 1 

I* + RI 

I + RI 

Cl 

2.0 x io -16 ( 4 . 2) - 1 

7 x io -17 (4.1) ±1 

I* + I 2 

I + I 2 

q 2 

1.9 x io* 1 1 ( 2 . 6) - 1 

1.9 x l O* 1 1 ( 2 . 6) - 1 

I* + I + RI 

I 2 + RI 

c i 

3.2 x io* 33 (3.2)* 1 

3.2 x io* 33 (3.2) ±1 

I* + I + I 2 

I 2 + I 2 

C 3 

8 x •) o -32 ( 1 . 8) ±1 

8 x io -32 ( 1 . 8) - 1 

R + RI 

r 2 + I* 

K 5 

3.2 x io -17 ( 3 . 2) - 1 

3.2 x io 17 ( 3 . 2 ) ± '' 

I* + RI 

RI 2 

k 6 

8.3 x 1 o* 1 8 ( 1 . 3 ) 15:1 

6.5 x io* 18 (1 . 1 ) ±1 


Reaction Rate for R + R ■* R 2 

The reaction rate coefficient for R + R -> R 2 is given for CF 3 radicals to be 

1.1 x IO * 11 cm 3 /sec by Hohla and Kompa (ref. 14). It is ejected to be slower for 
C 3 F 7 because of steric factors, although the additional degrees of freedom tend to 
stabilize this two-body reaction. Ibis value was likewise adopted by Turner and 
Rapagnani (ref. 15). Further studies were made by Kuznetsova and Maslov (ref. 16) on 
the normal and isomeric forms of the radical and with varying degrees of saturation 
of the gain medium. They find a value of 3.19 x 10 -12 cm 3 /sec within a factor of 

3.8 for the normal form and 7.7 x 10 -13 cm 3 /sec within a factor of 3.1 for the 
isomeric form. Later studies using transverse electron discharge excitation yielded 

4.2 x IO -12 cm 3 /sec for i-CjF^ (ref. 17). Values for recombination are given by Fisk 
(ref. 18) as 6 x 10 -13 cm 3 /sec for n-C 3 F 7 and 1.5 x io* 12 cm 3 /sec for i-C 3 F 7 , 
although his logic in arriving at these values is not detailed. The uncertainty in 
Kuznetsova and Maslov's original estimates was greatly reduced in their later work 
(ref. 19) with no substantial disagreement with their earlier values (ref. 16). 
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Later measurements of Skorobogatov and Slesar (ref. 20) yielded the value of 
1.77 x 1 0 -1 2 cm 3 /sec for i-C-jFy The values adopted for this study are shown in 
table VII. 


Reaction Rate for R + I -> Rl 

The principal reaction by which the lasant material is recovered is R + I -* RI. 
The rate proposed by Hohla and Kompa (ref. 14) and used by Turner and Rapagnani 
(ref. 15) is 5 x 1 0 -1 1 cm 3 /sec. The value obtained by Kuznetsova and Maslov 
(ref. 16) is 2.8 x lO - ^ 1 cm 3 /sec for both n-C^F^ and i-C^F^ within a factor of 3.1. 
Hie value found according to Skorobogatov et al. (ref. 21) is 1.5 x 10~ 10 cm 3 /sec 
within a factor of 3, which is equal to the value found by Beverly and Wong 
(ref. 17). The more recent values of Kuznetsova and Maslov (ref. 19) are somewhat 
lower than their earlier value. Their rate coefficients are 0.7 x 10 -11 cm 3 /sec for 
n-C^Fy and 0.6 x 10 cm J /sec for i-C 3 F 7 with an uncertainty factor of 1.15. A more 
recent value of Skorobogatov and coworkers (ref. 20) is 3.1 x 10 -11 cm 3 /sec for 
n-C 3 F 7* ' Ihe value 2,3 x 10-11 cm 3 /s ec is adopted for n-C 3 F 7 with an uncertainty 
factor of 3.5, and the value 3.9 x 10 -11 cm 3 /sec is adopted for i-CjFy with an 
uncertainty factor of 4.3, as shown in table VII. 


Reaction Rate for I + I ■* I 2 

Ihe reaction I + I -> I 2 is important because it competes with the regeneration 
of the lasant gas, and most importantly, the product is a serious quencher of the 
upper laser level (ref. 22). The iodine atoms approach along an attractive potential 
in the 1 E state but must be stabilized by collision with a third body. Ihe 3 it state 
has a weak repulsive potential allowing for close approach of the two iodine atoms 
with collisional transitions to the 1 Z state being likely through a third body. The 
efficiency of the stabilization is dependent on the mass ratio of the third body to 
the iodine atom and the van der Waals force of the third body. There is a rough 
inverse correlation between boiling point of the third-body material (that is, strong 
van der Waals potential) and the recombination coefficient (ref. 23). 

The most important third body for solar pumping is the alkyl iodide molecule, 
for which the rate coefficient is taken as 4.5 x 10 -J cm /sec by Turner and 
Rapagnani (ref. 15). Measurements of Kuznetsova and Maslov (ref. 16) yielded 
8.5 x 10- 32 cm 6 /sec for n-C 3 F 7 I and 8.0 x 1(T 32 cm 6 /sec for i-C-jFyl. The value of 
3.9 x 10 -2 ^ cm^/sec is given by Ershov et al. (ref. 24). 

The second most important third body in the gas is the I 2 molecule. The rate 
coefficient for this reaction is reasonably approximated as 3.8 x 1 0 -30 cm®/sec 
within a factor of 1.3 as established from references 15, 25, 26, and 27. However, 
this molecule must always give only a minor contribution because I 2 concentrations 
remain small for laser operation. 


Reaction Rate for R + RI ■* R 2 + I 

The reaction R + RI -> R 2 + I is hindered by steric factors and activation energy 
on the order of 2 eV. Measurements of Andreeva et al. in CF.,1 (ref. 28) yield a 

<1 fz o *5 

value less than 3 x 10 cm J /sec. 


15 



Reaction Rate for R + I* ■* RI 


The reaction R + I* •> RI is hindered by steric factors and a repulsive 

potential. Hie value measured for CF-. by Andreeva et al. (ref. 28) is 
*”123 ^ 

3.5 x 10 cm /sec, and a somewhat smaller value is expected for C 3 F 7 as a result 
of steric factors. Kuznetsova and Maslov (ref. 16) arrive at 1.12 x 10“^ 2 cm 3 /sec 
within a factor of 2.9 for n-C 3 F 7 and 0.5 x 10 ~ 12 cm 3 /sec within a factor of 7 for 
i-C ? F 7 . Hie transverse-electron-discharge-excited gas of Beverly and Wong (ref. 17) 
indicates a value of 2.7 x 10 -12 cm 3 /sec for i-C 3 F 7 « Later work of Kuznetsova and 
Maslov (ref. 19) yields the somewhat lower values of 9 x 10 -14 cm 3 /sec for n-C 3 F 7 and 
4 x 10 -14 cm 3 /sec for i-C 3 F 7 . More recently, the value of 6 x 10 “ 13 cm 3 /sec for 

i-C 3 F y is given by Skorobogatov and Slesar (ref. 20) , and a somewhat larger value of 

3 x 1 0 - ”* 2 cm 3 /sec is found by Ershov et al. (ref. 29). Vinokurov and Zalesskii 

(ref. 30) suggest the rate coefficient for i-C 3 F 7 to be 1 * 10 cm^/sec. Hie value 

of 5.6 x 1 0 _1 3 cm 3 /sec with an uncertainty factor of 6.2 is assumed for n-C 3 F 7 , and 
1.7 x 10 “ 13 cm 3 /sec with an uncertainty factor of 17 is assumed for i-C 3 F 7 , as shown 
in table VII. 


Reaction Rate for I* -> I 

Loss of electronic excitation energy through processes other than laser oscilla- 
tion is in direct competition with power generation. The loss of electronic energy 
can be through several processes. Thermal collisions are nearly adiabatic, so that 
energy transfer to translational motion is quite slow (ref. 31). The main modes of 
deexcitation will be energy transfer to vibrational-rotational levels and complex 
formation. 

The deexcitation by the lasant molecules proceeds through exciplex formation 



RI + I + hv 


+ I 


with subsequent decay through radiation or energy transfer within the complex. Hie 
dagger superscript denotes vibrational excitation. Hie branching ratios have been 
studied (refs. 13 and 32) by observing the exciplex emission. Hie overall rate for 


RI + I* -> RI + I 


has been measured by several groups. Turner and Rapagnani (ref. 15) recommend 
8 x 1 o -1 ® cm 3 /sec. Somewhat smaller values were found by Kuznetsova and Maslov 


(ref. 16) as 2.3 x 10 - ^® cm 3 /sec for n-C 3 F 7 I and 2.8 x 10 -1 ® cm 3 /sec for i-C 3 F 7 I. 
Accurate studies of chemiluminescence b^ Stephan and Comes (ref. 33) and Comes 


and Pionteck (ref. 34) yield 4.6 
2.0 x 10~ 16 cm 3 
value, 1.7 x 
shown in table VII. 


x 10 


cm /sec for n-C 3 F 7 I. 
m J /sec, was found by Dobychin et al. for n-C 3 F 7 I, 
10 -1 ^ cm 3 /sec, was found for i-C 3 F 7 I (ref. 35). 


A larger value, 
and a much lower 
Hie values adopted are 
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The deexcitation by I 2 is probably through the formation of a triatomic complex 
(ref. 36). The value adopted by Turner and Rapagnani (ref. 15) was probably influ- 
enced by the early value (5 x 10 -12 cm^/sec) given by Donovan and Husain (ref. 37). 

A larger value of 3.6 x 10 -11 cm 2 /sec is given by Beverly and Wong (ref. 17). 

Kartazaev et al. (ref. 38) find 2.1 x lO - ^ 1 cm 2 /sec. Most of the more recent values 

approach 3 x 10 -11 cm 2 /sec (refs. 27, 36, and 39). 


Reaction Rate for RI + RI •> R 2 + I 2 

Although the enthalpy change is easy to achieve, RI + RI 4 + I 2 is an unim- 

portant reaction, inhibited by both steric factors and a high activation energy. The 
rate is assumed to be zero. 


Reaction Rate for I* + I I, 


The reaction I* + I •> X 2 is important because it is the main source of I 2 gen- 
eration before the medium is saturated. The 3 it state is attractive and easily recom- 
bines with a third body present to absorb relative translational energy in forming 
the bound state. The most important third body is the lasant gas, for which coef- 
ficients have been measured. Turner and Rapagnani (ref. 15) give a value of 

x 1( T 32 cm®/sec. More 

shown in 


.-33 


1 x icT~’~' cm 6 /sec, and Beverly and Wong (ref. 17) use 1.56 x 10 -32 cm 6 /sec. 
recently, Stephan and Comes (ref. 33) give the value 1 x 10 -32 cm 6 /sec. As 
table VII, the value assumed here is the geometric midpoint 3. 2 x 10 -33 cm^/sec with 
an uncertainty factor of 3.2, which is sufficient to span the experimental values. 


The reaction as stabilized by the I 2 molecule is given as 4.3 x 10 32 cm 6 /sec 
by Turner and Rapagnani (ref. 15) and 1.48 x 10 -2 ^ cm^/sec by Beverly and Wong 
(ref. 17). The value of 8 x 10 -22 cm®/sec with an uncertainty factor of 1.8 is 
assumed in table VII. 


Reaction Rate for R + RI -*■ R 2 + I* 

The reaction R + RI -* R 2 + I* was proposed by Andreeva et al. (ref. 40) as a 
postphotolysis metastable production mechanism. The reaction is hindered by steric 
factors and a high activation energy and is most likely to proceed through hot 
radical reaction. Previous studies of Andreeva et al. (ref. 28) indicate this reac- 
tion to be unimportant. Direct measurements of Zalesskii and Krupenikova (ref. 41) 
yield the rate constant 3.2 x 10 - ^ cm 3 /sec within a factor of 3.2. This reaction 
was proposed by Hohla and Kompa as a possible explanation for the postphotolysis gain 
increase observed in their experiments. 


Reaction Rate for I* + RI RI* 

The RI and I* system forms a weakly bound state (-0.07 eV for n-CjP^I and 
-0.04 eV for i-C 3 F 7 I) through three-body collision (refs. 13 and 32). The exciplex, 
once formed, decays through emission of 1.3 pm radiation or through transfer of the 
electronic excitation energy into vibrational-rotational motion of the RI molecule. 
The overall effect is loss of atomic iodine metastables or quenching of the 
metastable state. The rate for passing through the exciplex radiative channel is 
measured by Gerck (ref. 32) to be 6.6 x 1 0 -1 8 cm 3 /sec for n-C 3 F 7 I and 
5.8 x 10 - 1® cm 2 /sec for i-C^Fyl. Somewhat higher values are observed by Smedley and 
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Leone (ref. 13) as 10.4 x 
1-C 3 F 7 I * 


10 -18 cm 3 /sec for n-C 3 F 7 l and 7.2 x io 


— 1ft ft 

cm J /sec for 


Other Reaction Rates 


The remaining reactions in 
no reaction rate constants were 


tables V and VI are not ejected to contribute, and 
found in the literature. 


SPECTRUM AND GAIN COEFFICIENT OF METASTABLE IODINE 

The photodissociation and subsequent chemical reactions govern the production 
and loss of inversion in the gas. The ability of the gas to amplify a light signal 
depends on the inversion and stimulated emission cross section. Hie spectral proper- 
ties and their relation to the stimulated emission cross section are discussed in 
this section. 

The valance shell of the iodine atom is a (5p) 3 configuration with total 
L = 1, S = 1/2, and a multiplicity of 2. The ground state is J = 3/2 with the 
first excited state due to spin-orbit interaction J = 1/2. In the Russell-Saunders 
scheme, Garstang (ref. 42) finds a magnetic transition rate of 7.8 sec -1 and a small 
electric quadrupole rate of 0.055 sec -1 which yield 0.127 sec as a lifetime of the 
^ P 1/2 state ' an< 3 transition energy of 7605 cm -1 . More recent calculations of O'Brien 
and / Bowen (ref. 43) yield 0.11 sec. Experiments were performed by Dervpnt and Thrush 
(ref. 44) by which the lifetime was measured relative to the 0 3 ( A) lifetime. 

Assuming the O^A) lifetime to be 3880 sec yields a lifetime of 0.17 + 0.04 sec for 

2p 

1 / 2 * 

1 ft 7 

Natural iodine occurs as 100-percent I , for which the nuclear spin is 5/2. 

Hie resultant total angular momentum of the whole atom is 


F = 2, 3 


(1 7) 


for the 



upper 


level and 


F = 1,2, 3, 4 


(18) 


for the 2 P;3/2 g roun( 3 state. Because of interaction between the magnetic field of the 
atomic orbital and the nuclear magnetic moment, the sublevels of different F values 
are split into distinct lines. In accordance with magnetic selection rules 
AF = 0, ±1 , there are six lines in the 2p 1 / 2 * 2p 3/2 transition. Hie energy split- 
ting (AE) of the upper levels calculated by Zuev et al. (ref. 45) is 


AE = -0.371 cm 
F=2 


AE„ = 0.265 cm 
F=3 
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and the splitting of the lower level is 


AE = -0.1 1 8 cm 

r = i 


AE = -0.094 cm 
F=2 


-1 


( 20 ) 


AE = -0.028 cm 
F=3 


-1 


AE = 0.113 cm 
F=4 


-1 


The energy differences between the levels are shown in table VIII along with results 
of experimental measurements (refs. 45 and 46). 


TABLE VIII.- HYPERFINE SPLITTING OF IODINE 2 P STATES 





AE„ „ 
F-| F 

, cm 1 

2 



F 1 - F 2 


2 P 

p 1/2 



2p 3/2 



Theory 

Experiment 

Ref. 

Theory 

Experiment 

Re f . 

1-2 




0.025 

0.026 ± 0.003 

45 






0.025 

46 

2-3 

0.636 

0.662 ± 0.010 

45 

0.0656 

0.068 ± 0.002 

45 



0.69 ± 0.02 

46 


0.065 

46 

3-4 




0.141 

0.141 ± 0.002 

45 






0. 1 42 

46 


Zuev et al. (ref. 45) calculated the transition rates for the six lines and 
obtained values shown in table IX. Ihe total lifetime of the upper state was found 
to be 0.13 sec, which is in good agreement with Garstang (ref. 42). 


The hyperfine levels are mixed by thermal collisions in a time which is short 
compared with its radiative lifetime of 170 msec. Ihe level densities are related by 
the equilibrium requirement 




( 21 ) 


19 


TABLE IX 


HYPERFINE TRANSITION RATES AND LINE STRENGTHS 


Transition 

Rate, 

-1 

sec 

3-4 

5.0 

3-3 

2.1 

3-2 

.6 

2-3 

2.4 

2-2 

3.0 

2-1 

2.3 


Relative intensity, percent 


Calculated 
(this study) 

37.9 

16.0 

4.5 

13.0 

16.2 

12.4 


Experiment 
(ref. 45) 

36.5 ± 1.5 
16.3 + 0.8 
6.2 ± 1.1 
13.0 + 0.6 

15.2 ± 0.7 

12.3 + 0.5 


Calculated 
(ref. 45) 

37.5 
16.2 

4.6 
13.0 
1 6.2 

12.5 


where gp and gp i are the degeneracies, and Ep and Ep, are the corresponding 
energy levels. In equilibrium, the forward and backward reaction rates are likewise 
equal 


K [l*l = K [i* 1 (22) 

f L F J b L F' J 


where and Kj-, are the forward and backward rate coefficients, respectively, so 

that 


K 


b 



(23) 


Since the energy separation between the hyperfine sublevels is small compared with 
the thermal energy, the temperature -dependent factor may be ignored. The ratio of 
level densities is assumed to be proportional to their degeneracy, the line intensi- 
ties are calculated, and the corresponding relative intensities for the transition 
rates are compared with experiments in table IX. 


Broadening of Spectral Lines 

There are several contributions to the width of a given atomic spectral line. 
First is the finite lifetime of the level resulting in energy uncertainty at emis- 
sion. This spectral broadening is impossible to observe under ordinary conditions. 
There is a Doppler shift due to the relative motion of the atom to the observer. For 
a gas in thermal equilibrium, the Doppler width is 

/ \ 1/2 

D /l-mo T — o \ ' 

V 

o 


v (|t/m 


Av d 


z ki v in x 


m 


Av q = 7.162 x 10 


-7 


20 


(24) 



where T is the temperature, and M is the mass of the gas molecule. The relation 
between the peak of the line shape and the line width is 


g o (°) 



(25) 


and applies at low gas pressures (ref. 47). 

At higher pressures, the interruption of the emission process by binary colli- 
sions causes random phase changes in the oscillator frequencies, which further 
broaden the emission spectrum (ref. 48) . Hie associated line shape is Lorentzian, 
and the line width is related to the gas collision frequency. Similar to collision 
broadening is van der Waals broadening, associated with the formation of weakly bound 
van der Waals molecules. The transition frequency then corresponds to the separation 
of the two potential curves at the relative nuclear separation at which the transi- 
tion occurred in accordance with the Franck-Condon principle (ref. 49). The colli- 
sion and van der Waals broadening are pressure dependent and approximately Lorentzian 
in profile. The line width is assumed linear in pressure as 

Av = a p (26) 

1 

where is the pressure-broadening coefficient. The peak of the line shape is 

given by 


g i (0) 


_2 

n Av 


(27) 


The line shape function exhibiting the combined effects of the Doppler and Lorentz 
components is the Voigt function. In the present work, the line shape is approxi- 
mated by a Lorentzian profile 


g( v) 


Av 

2n 


1 

(v - v ) 2 + (Av/2) 2 
o 


(28) 


since overlap of the pure Doppler profile of adjacent lines is negligible. The line 
width is taken as 


Av = a + a p (29) 

o 1 

where p is the alkyl iodide pressure, and a Q is related to the Doppler line width 
in equation (24). 
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Gain and Stimulated Emission Cross Section 


The inversion density on a given line is expressed as 


AI 


FF' 



( 30 ) 


for which the corresponding stimulated emission cross section is 


X 2 a 


FF 


. ( V) = 


FF' 


8n 


'FF 


,(v) 


(31) 


where g , ( v) is the profile function of the F > F' transition, and A p p, i s the 
transition rate. The gain coefficient at a frequency v is 


G(V) a FF' (V) AI FF- (32) 

FF' 

which has been studied by several groups. Zuev et al. (ref. 45) measured the line 
width of the iodine luminescence directly and measured the gain on the 3-4 transi- 
tion. They conclude that the A 34 transition rate is 3.4 ± 1.1 sec -1 / which implies 
an Einstein coefficient of 


A = 5.4 ± 2 sec -1 (33) 

in good agreement with the value 6.22 ± 1 .46 sec -1 of Derwent and Thrush (ref. 44) 
but somewhat lower than the 7.86 sec -1 of Garstang (ref. 42). The stimulated emis- 
sion cross section on the 3-4 transition was likewise measured by Fuss and Hohla 
(ref. 50) . The line width was obtained from the measured stimulated emission cross 
section as 


Av 


A 34 X 


a 4iz 

34 


(34) 


— 1 

where they assumed the value of A 34 to be 5.0 sec according to the calculations 

of Zuev et al. (ref. 45). They obtain for a of equation (29) the value of 

20.2 MHz/torr. Baker and King (ref. 51) measure a 34 at 100 torr to be 

0.9 ± 0.2 x 10 -18 cm 2 , which is in good agreement with Zuev et al. , who assumed 

that A^^ = 3.4 ± 1.1 sec -1 . Nearly perfect agreement between Zuev et al. and Baker 

and King is obtained for a 34 if A 34 = 3.7 sec -1 is used with the Zuev et al. 

values. 
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A detailed study was undertaken by Mukhtar, Baker, and King (ref. 52) in which 
line width was measured independently of the cross section, and the line- 

broadening parameter a was found to be 15 ± 4 MHz/torr. Results obtained by 
Babkin et al. (ref. 53) are quite distinct from the results of references 44, 50, 51, 
and 52, although the Babkin et al. broadening parameter for CF^I is in line with Fuss 
and Hohla (ref. 50) . The pressure-broadening parameters of various groups are listed 
in table X along with estimated temperatures at the time of measurement. The room- 
temperature parameter is estimated assuming a T dependence. The value adopted here 
is 14.8 ± 4 MHz/torr. 


TABLE X.- PRESSURE-BROADENING PARAMETER FOR C 3 F ? I GAS 


a-, (T), 

T, 

a! (298 K), 

Reference 

MHz/torr 

K 

MHz/torr 


18 + 4 

530 

13.5 ± 3 

Zuev et al. (45) 

20 

420 

16.9 

Fuss and Hohla (50) 

1 5 ± 4 

350 

13.9 ± 3.7 

Mukhtar et al. (5 2) 

Adopted value 


14.8 ± 4 



EFFECTS OF DIFFUSION 

Diffusion of the iodine metastable within the laser tube with finite dimension 
is an additional process by which inversion is lost within the gas. The diffusion 
proceeds in two steps associated with therma lization of the hot photodissociation 
products followed by slow thermal diffusion. 

The iodine metastable is produced with kinetic energy of 0.79 eV and undergoes 
4 to 5 collisions before thermalization takes place. The distance traveled in ther- 
malization is 


d ~ 0.5 mm-torr 


(35) 


which is small compared with the radius of the laser tube (»3.5 mm) even at 1 torr. 
However, the thermal diffusion may be quite important. Considering the expected high 
sticking factor of the iodine atom at the quartz wall (ref. 54), such collisions 
appear as a loss of excitation energy from the inversion density. 

The particles produced at an instant in time at the origin of a coordinate frame 
as a result of gaseous diffusion are distributed in space at a later time t accord- 
ing to 


n 


I* 


(r, t) 


2n 

o 

7t( 4Dt) 3/2 



(36) 
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where n Q is the initial density of metastables, r is the distance from the ori- 
gin, and D is the diffusion coefficient (ref. 55). In the laser, the source of I* 
is relatively uniform over several millimeters, so that losses due to diffusion from 
one point are compensated by diffusion from adjacent source points. Only by diffu- 
sion to the glass surface are particles lost. Density changes on the tube axis are 
then proportional to 


n x *(0,t) = 


„ /•“> r r 

2n I o 

-r " 

J 


r dr 


1 


(4Dt) 


3/2 


exp 


2 2 
- (r + JT) 


4Dt 


n (0,t) = n 1 
I* o 



(37) 


where r Q is the inside radius of the laser tube. A diffusion time constant can be 
defined by using equation (37) as 


x 


D 



(38) 


where x is the mean diffusion time from the centerline to the tube wall. 

D 

The diffusion coefficient for isotropic scattering is related to the elastic 
scattering mean free path \ as 



(39) 


where v is the mean thermal speed of the diffusing particle. The iodine atomic 
radius is 215 pm (ref. 9), the alkyl iodide radius is estimated from binary diffusion 
(ref. 56) to be 400 pm, and the mean thermal velocity of the iodine atom is 
1.92 x 1 o 4 cm/sec (298 K). Using the iodine atomic radius and the alkyl iodide 
radius yields 3.8 x 10 -1 ^ cm 2 for the scattering cross section and 0.0076 cm-torr for 
mean free path so that equation (39) yields 


2 

D ~ 73 cm -torr/sec 


(40) 



The laser tube radius of 3.5 mm yields a diffusion lifetime of 


8.5 


-4 

x 10 sec/torr 


( 41 ) 


for the present experiments. 


PHOTOLYSIS OF ALKYL IODIDES 

The chemical and physical processes being outlined in previous sections, a math- 
ematical model describing the rate of change of specific densities of constituents in 
the laser tube is described in this section. The major kinetic pathways are first 
summarized. The initial excitation of the gas is through photodissociation into 
excited iodine and the C 3 F 7 radical denoted by the zigzag line in figure 7. The loss 
of excited iodine through radiative decay is slow (see zigzag line between I* and I 
in fig. 7), so that excited iodine loss is mainly through the slow recombination 


R + I* -> RI 


which allows large metastable populations. Further losses of excited iodine occur 
through collisional quenching, especially if contamination of the gas with I 2 , H 2 0 ' 
or C >2 occurs. The ground-state iodine rapidly combines with the free radical C^Fy. 
This recombination is necessary for maintaining the inversion. More slowly, the 
ground-state iodine recombines with either ground-state iodine or metastable iodine 
to form iodine molecules. Due to the presence of intense visible radiation, the 
iodine molecule photodissociates to form ground-state and metastable iodine. 




Figure 7.- Major kinetic pathways of the iodine laser. 
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The major kinetic pathways are summarized in figure 7. The photodissociation 
rates are given by equation (15), where the parameters are those in table II. The 
maximum gain will occur on the tube centerline on which the kinetic equations are to 
be solved. Diffusion effects are ignored except for diffusion of the iodine atoms to 
the inner wall, where they are assumed to have a sticking factor of unity (ref. 54) . 

The RI molecule photodissociates at a rate ^ to form R and I*. The hot 
radical R is assumed to collide with the RI molecule to form dimer R 2 and ground- 
state iodine with rate coefficient K 4 » The RI is reformed by recombination with I or 
I* as depicted in 


d [RI] 
dt 


= K 1 [R] [I* ] + K 2 [R][I] - K 4 [R][RI] 




(42) 


where brackets around a chemical species denote its density. Note, this is an 
approximate representation for hot radical reaction. The source of the free radical 
R is photodissociation of RI. Losses occur through recombination of R with R, I, or 
I*. Additional loss occurs through hot radical reaction 


d [R] = £ 

dt S 


- [R] [I * ] - K 2 [R] [I] - 2K 3 [R] 2 


- K [RI] [R] 
4 


(43) 


The dimer is formed by radical recombination or hot radical reaction 


d[R 2 ] 

dt 


= K 3 [R ] 2 


+ K 4 [R] [RI] 


(44) 


The iodine molecule is formed through three-body recombination of iodine atoms 
and metastables. The metastables do not recombine because of mutual repulsion. The 
iodine molecules are lost through photodissociation. The rate of change of iodine 
molecules is given by 


d t ! 2 1 2 2 

= C 1 [I] [I*] [RI] + C 2 [I] [RI] + C 3 [I] [I*] [I 2 ] + C 4 [I] [I 2 ] - ? 2 (45) 

The iodine metastable is formed through photodissociation of RI and 12* Ik is lost 
mainly through recombination with R. Slower loss processes include three-body recom- 
bination with iodine atoms and quenching with RI and I 2 molecules. Further loss 
occurs through radiative decay and diffusion to the cell wall. The resultant rate 
equation for metastable density is 


-5 + l 2 ~ K [R][I*] - C [I] [1*3 [RI] - C [I] [I*] [I ] 

[I*] 

- Q 1 [I*][RI] - Q 2 [I*] [ I 2 ] - A [ I* ] 


(46) 
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Ground-state iodine atoms are produced through photodissociation of I 2 , quenching, 
and radiative decay of iodine metastable. A small source is hot radical reaction 
with RI. Losses are mainly through R + I recombination, and smaller losses appear 
through three-body recombination and diffusion to the cell wall. The total atomic 
iodine change in time is given by 


d [ I ] 
dt 


= 5 2 [I 2 1 + Q [I*][RI] + Q 2 [I*][I 2 ] + K 4 [R][RI] 


+ A [ I* ] - C [I][I*][RI] - 2C 2 [I] 2 [RI] - C 3 [I] [I*] [ 1 2 3 

— 2C 4 [I] 2 [X 2 ] - K 2 [R][I] - (47) 


Equations (42) through (47) apply at any location on the laser centerline if no lon- 
gitudinal diffusion is assumed. Kinetic rate coefficients are given in table VII. 

Consider the initial illumination of the lasant gas, for which all other con- 
stituent densities are initially zero. The most important processes are 


RI + hv a- R + I* 


R + I* -> RI 


and to a lesser extent 


I* + RI -> I + RI 


and diffusion at low pressures. To a good approximation, the R and I* densities are 
equal, so equation (46) becomes 


d ^ - ] = ^ - K [I*] 2 - [RI] [I*] - (48) 

' D 


the solution of which is given by 


[I*] 


-2 

- T 


2K 


1 


-1 

A + X 


1 - exp (-At) 

+ ( A - t 1 ) exp (-At) 


(49) 


27 



where 


A 



-2 

+ T 


(50) 


and 


t _1 = Q 1 [RI] + t" 1 


(51) 


It is clear that the inversion density is closely related to the rate coefficients 
Ki and Q-] and to the diffusion constant Td . 

The time constant associated with reaching steady state is found from equa- 
tion (50) to be 


A 1 - ^ (msec) 
VPC 


(52) 


where C is the solar concentration, and p is the alkyl iodide pressure in torr. 

In actual practice, the light input is modulated by a chopper for which the rise time 
is on the order of a few milliseconds. (See fig. 8. ) Ihe effective time constant is 
then extended from the value in equation (52) by the amount 


-1 100 
A ~ — 

^pC f (t) 


(53) 


where f(t) is the light modulation function, which is initially zero and which 
rises to unity in «2.5 msec. Ihe solution is more accurately given by 


[I*] 


A - t~ ^ J - exp ( -At/3 ) 

2K^ 1 + exp(-3At) 


(54) 


where A 


is 


time dependent according to equation (53) 




Time, msec 


Figure 8.- Light intensity modulation function used in model compared 
with light measurements for chopper used in experiments. 


In most cases of interest, threshold is reached in a fraction of a millisecond, 
for which equation (54) is far from its steady-state value, so that expansion of 
equation (54) for small values of At yields 


[I*] 




(55) 


as one would expect on observing equation (48) . In equation (55) , t is the rise 
time of the chopper modulation function (t c ~ 2 msec). Hie relation of the I* den- 
sity to the pump source yields the spatial dependence of the inversion density for 
which the average gain over the tube length may be found. 


MAXIMUM GAIN 

The equations of photolysis are integrated numerically at the most intense posi- 
tion of the light source. The pump light is modulated by the mechanical chopper 
according to a modulation function which rises in 2.5 msec, holds near unity for 
4 msec, and declines over the next 2.5 msec. The maximum gain occurs on the 
3-4 transition because of the high population in the F = 3 upper level (eq. (21)) 
and the high transition rate to the F = 4 lower level (table IX). In accordance 
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with the definitions of Fuss and Hohla (ref. 50) , the gain coefficient for the 3-4 
transition is taken as 


G 


34 


Z a (v ) Al 

FF' 34 FF' 

FF' 


G 


34 


^ > °FF' ^ V 34^ 
FF' 



(56) 


where the total population inversion between P 1 and p 3/2 


Al = [I*] -If!- (57 

and the ratio 83/^2 + 93 ) is the approximate fraction of the inversion which con- 
tributes to the 3-4 transition (ref. 50). The stimulated emission cross section is 
then taken as 


0 34 



FF 


(V 34 ) 


(58) 


The individual values of ° FF i( v 34 ) are evaluated using equation (31). The inver- 
sion varies along the axis of the tube in approximate proportion to the illumination 
by the arc image. (See eqs. (4) and (55).) The total round-trip amplification of 
the gas on the axis of the tube is then 
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which must balance the losses at the tube windows and end mirrors (with effective 
reflectances R-| and R 2 ^ as 


R 1 R 2 


exp 2 


T7fi Lo 3*4 ai 


(60) 
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The inversion density Al depends mainly on the kinetic coefficients and 

and the diffusion coefficient for I* in RI as seen in equations (49) to (51). The 
stimulated emission cross section depends on the pressure -broadening coefficient 
a . These parameters will now be investigated from the point of view of threshold 
measurements on solar-simulator-pumped lasers. 


COMPARISON WITH EXPERIMENTAL RESULTS 

Threshold was investigated experimentally with the setup depicted in figure 1. 
Before each experimental run, the laser tube was evacuated to “10 - ® torr and filled 
with purified lasant gas (n-C^F^l) through a copper gas-handling system. Prior to 
use in experiments, the lasant was purified by distillation. Because of the temporal 
oscillations of the arc light source, many experimental runs were generally made at 
each fill pressure in an attempt to identify ideal exposure conditions. The identi- 
fication of maximum pump rate was made by placing threshold measurements at each 
pressure into a sequence of decreasing values, which should provide a sequence which 
approaches the model from above. Threshold was identified by the beginning of the 
first laser pulse and recorded as the fraction of the peak light intensity of the 
solar-simulator output pulse. Note that the threshold condition expressed in this 
way depends on the chopper geometry and speed, since threshold is achieved before the 
photolysis establishes equilibrium with the solar-simulator light source. Before an 
attempt to compare the model with the experiments was made, the sensitivity of 
threshold to the coefficients , q 1 , and aq was determined. First, the kinetic 
equations were integrated numerically by using the mean values of all coefficients 
shown in table VII. The threshold as a function of pressure for the mean values of 
parameters is shown in figure 9 as the lower curve. Then each of the 



Figure 9.- Mean threshold and threshold sensitivity to experimental uncertainty 
in specific rate coefficients and pressure broadening. 
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parameters f Q-j , and cc-\ was allowed to take its maximum value specified in 
table VII one at a time to determine the sensitivity of threshold to each. These 
results are shown in figure 9. It is seen that affects the threshold almost 

equally at all pressures. causes little effect below 5 torr but a very large 

effect at 60 torr. Hie pressure-broadening coefficient a 1 shows a slight effect at 
low pressure and greater effect with increasing pressure. Parameters were so chosen 
as to provide the most reasonable agreement with experiments from the three conic 
collectors. Results are shown in figures 10, 11, and 12. Hie calibration shown in 
figure 2 for the cone with a 2 3.6-cm-diameter aperture results in CqIIl = 618, which 
compares with a more recent calibration on July 31, 1981, which yielded » 540. 

The results shown in figure 12 are calculated for the July 31, 1981, calibration. 

Hie final values used were 


-12 3 

= 2.0 x io cm /sec 
-16 3 

Q = 8.4 x 10 cm /sec 
=17.6 MHz/torr 


It is noted that Q 1 and a 1 are similar in their effect on threshold. Hie high 
values of Q-, used here may be caused by gas impurity. These results are summarized 
in table XI along with results of the previous analysis of the solar-simulator-pumped 
iodine laser experiments (ref. 4). 



Figure 10.- Comparison between theoretical and e>qperimental threshold for 

cone with 9 .6 -cm -diameter aperture. 
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TABLE XI.- IODINE LASER KINETIC RATE COEFFICIENTS ACCORDING TO 
ANALYSIS OF SOLAR-SIMULATOR EXPERIMENTS WITH n-C 3 F ? I 

The factor in parentheses gives the uncertainty limits 
associated with the coefficient 


Coefficient 

From table VII or X 

Solar simulator 

K 1 , cm 3 /sec 

5.6 x 10“ 13 (S^)* 1 

2.0 x 1 0 -1 2 

K 2 / cm 3 /sec 

2.3 x 10 -1 1 (3.5)* 1 

4.3 x 10' 11 (1 .O)* 1 

K^, cm 3 /sec 

2.6 x 10 -12 (4) ±1 


K 4 , cm 3 /sec 

3 x 1 o -1 6 


Q 1 , cm 3 /sec 

2.0 x 10 -16 (4.2) ±1 

8.4 x 10" 16 

q 2 , cm 3 /sec 

1.9 x 1 0 _1 1 ( 2 . 6 ) — 1 

3.1 x 1 0 -1 1 < 1 . 6) * 1 

C 1 , cm 6 /sec 

3.2 x IQ' 33 (3.2)* 1 


C 2 , cm 6 /sec 

8.5 x IQ' 32 (5.3)* 1 

3.7 x icT 32 (2.3) ±1 

Cj» cm^/sec 

+1 

CO 

• 

t— 

m 

l 

o 

X 

00 


C 4 , cm 6 /sec 

w 

• 

00 

X 

o 

! 

u> 

o 

• 

u> 

H- 


, MHz /torr 

14.8 ± 4 

17.6 


CONCLUDING REMARKS 

A detailed model of the chemical kinetics of the n-C-^Fyl solar-simulator-pumped 
iodine laser is utilized to study the major kinetic processes associated with the 
threshold behavior of this experimental system. Excited-state diffusion to the cell 
wall is the dominant limiting factor below 5 torr. Excited-state recombination with 
the alkyl radical (R) and quenching by the parent gas control threshold at higher 
pressures. Treatment of the hyperfine splitting and uncertainty in the pressure 
broadening are important factors in fixing the threshold level. 

The trends of the experimental threshold measurements are well represented by 
the present model. Even the scaling with respect to peak light input intensity and 
gain length for the three conic collectors is quite consistent. The final value of 
the recombination rate for 


I* + R -> RI 


is on the high side of the results reported by most other workers. Compared with 
values for i-C 3 F 7 recombination, our value is slightly smaller than the recent value 
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of Ershov et al. (Soviet J. Qaantum Electron., vol. 8, no. 4, 1978) but 3 times 
larger than Skorobogatov and Slesar (Vestn. Leningr. Univ. Fiz . Khim., vol. 4, no. 1, 
1979). lhe metastable quenching rate is quite high and may be indicative of contami- 
nation of the parent gas. lhe final value of the pressure -broadening coefficient is 
large compared with others. There is some dependence of pressure broadening and 
quenching coefficient so they can, to some degree, be traded off without affecting 
threshold in the model. If the gas is sufficiently contaminated, then pressure 
broadening may be less than indicated. 

The present compilation of data indicates that the isopropyl form of the alkyl 
iodide should have a lower threshold because of the larger photodissociation rates 
and slower recombination and quenching rates. It would be interesting to make such a 
comparison. Finally, experiments in which better control over experimental light 
input and gas contamination effects would be useful in resolving disagreement among 
various experimenters on the recombination and quenching rates. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 12, 1983 
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